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Abstract 

Transcription factors are proposed as suitable targets for the control of traits such as yield or food quality in plants. 
This study reports the results of a functional genomics research effort that identified ATHB17, a transcription factor 
from the homeodomain-leucine zipper class II family, as a novel target for the enhancement of photosynthetic capac- 
ity. It was shown that ATHB17 is expressed natively in the root quiescent centre (QC) from Arabidopsis embryos and 
seedlings. Analysis of the functional composition of genes differentially expressed in the QC from a knockout mutant 
(athb17-1) compared with its wild-type sibling revealed the over-representation of genes involved in auxin stimulus, 
embryo development, axis polarity specification, and plastid-related processes. While no other phenotypes were 
observed in athb17-1 plants, overexpression of ATHB17 produced a number of phenotypes in Arabidopsis including 
enhanced chlorophyll content. Image analysis of isolated mesophyll cells of 35S: - ATHB17 lines revealed an increase 
in the number of chloroplasts per unit cell size, which is probably due to an increase in the number of proplastids 
per meristematic cell. Leaf physiological measurements provided evidence of improved photosynthetic capacity in 
35S: - ATHB17 lines on a per unit leaf area basis. Estimates of the capacity for ribulose-1 ,5-bisphosphate-saturated and 
-limited photosynthesis were significantly higher in 35S: - ATHB17 lines. 

Key words: Chloroplast number, chloroplast size, homeodomain-leucine zipper, leaf development, photosynthetic capacity, 
transcription factor. 



Introduction 

Homeodomain-leucine zipper (HD-Zip) proteins fall into four coupled to a leucine zipper motif, which is thought to facili- 
classes (I-IV), and all contain a DNA-binding homeodomain tate homo- or heterodimerization (Ariel et al, 2007; Ciarbelli 
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et al, 2008; Elhiti and Stasolla, 2009; Harris et al, 2011). 
HD-Zip class II proteins are implicated in mediating the 
response of plants to a range of environmental conditions 
as well as in regulating developmental processes (Ariel et al. , 
2007; Ciarbelli et al, 2008; Elhiti and Stasolla, 2009; Harris 
et al, 2011; Ruberti et al, 2012). Five HD-Zip class II genes 
are known to respond to changes in light quality (Ciarbelli 
et al, 2008). For the remaining class II HD-Zip genes, lit- 
tle is known about their regulation other than that they are 
not responsive to changes in light quality (Ciarbelli et al, 
2008). The light responsiveness of members of the y subclass 
(ATHB2, ATHAT1, and ATHAT2) has been investigated the 
most thoroughly. For example, when grown in white light, 
ATHB2 expression is induced by shading (Steindler et al, 
1997). ATHB2 and ATHAT1 may play roles in the shade- 
avoidance response, as their expression is induced by treat- 
ment with far-red-rich light, implicating a possible role for 
phytochrome in the regulation of these two genes (Ciarbelli 
et al, 2008). 

Light quality and intensity can also affect plastid num- 
ber. For example, in plants grown under high-intensity white 
light, chloroplast number increases with a relative decrease 
in plastid size (Butterfass, 1979). Chloroplast division is a 
complex process that involves multiple distinct steps, which 
are controlled by different nuclear genes (reviewed by Pyke, 
1999; Osteryoung and McAndrew, 2001; Aldridge et al, 
2005; Pogson and Albrecht, 2011). Higher-plant chloroplasts 
develop in young mesophyll cells from small plastid initials 
called proplastids. Proplastids are present in meristematic 
cells and divide in association with cell division. This division 
ensures that a constant proplastid number (about 14; Pyke, 
1999) is maintained in all meristematic cells. As post-meris- 
tematic cells differentiate and expand, proplastids increase in 
size, develop into chloroplasts, and undergo approximately 
three rounds of division to produce the final complement of 
chloroplasts in a fully expanded mesophyll cell (Pyke, 1999; 
Aldridge et al, 2005). 

Here, we present results obtained with ATHB17 
(AT2G01430), which encodes a member of the a subclass 
of the class II HD-Zip protein family of transcription fac- 
tors (Riechmann et al, 2000; Ariel et al, 2007; Ciarbelli 
et al, 2008; Elhiti and Stasolla, 2009; Harris et al, 2011). 
This study describes the expression pattern of ATHB1 7 and 
characterizes the phenotypes of ^47W,B77-ovei"expressing 
plants and athbl7 knockout plants. The results describe 
the application of ATHB17 to regulate leaf morphology 
and chlorophyll content, the number and size of chloro- 
plasts, and the capacity for photosynthesis in Arabidopsis 
leaves. 



Materials and methods 

Plant material and growth conditions 

ATHB1 7-overexpressing lines were created in the Arabidopsis 
thaliana Heynh. ecotype Col-0. The ATHB17 (AT2G01430) clone 
was PCR amplified from a diversified cDNA library created using 
mRNA from different Arabidopsis tissues. This clone contained 
an isoleucine substituted for methionine at position 227 and also 



contained 64 bp of the native 5'-untranslated region (5'-UTR) 
immediately upstream of the start ATG and 116 bp of the native 
3'-UTR immediately downstream of the stop codon. Arabidopsis 
plants were transformed by the floral dip method (Bechtold et al. , 
2003) using Agrobacterium carrying a transformation construct 
containing a kanamycin resistance gene driven by the nopaline 
synthase (NOS) promoter and the ATHB17 clone downstream 
from the cauliflower mosaic virus 35S promoter. A control line 
generated by transforming Col-0 with an empty transforma- 
tion construct was used to determine the effect of ATHB17 
overexpression. 

Between 20 and 40 independent primary transformants were iso- 
lated on selection medium and transplanted into soil. For bulk seed 
production, plants were grown under 24 h light, in growth rooms 
(80-110 umol nr 2 s" 1 , 22 °C). Transformants were PCR geno- 
typed to confirm that they harboured the correct transgene using 
forward (5'-GCAAGTGGATTGATGTGATATC-30 and reverse 
(S'-TGAATTCTTCATCGTCTCCGTCTTC-S') primers. The 
majority of primary transformants showed alterations in leaf shape 
and dark green coloration, and in some instances reductions in plant 
size and alterations in inflorescence morphology were apparent. 
Expression of A THB1 7 mRNA transcripts was verified by reverse 
transcription-PCR (RT-PCR) on RNA extracted from leaves of 
40-d-old plants (see below). Lines that showed overexpression of 
ATHB1 7 by RT-PCR versus the control line were selected and used 
in subsequent experiments. 

For morphological analysis and photosynthetic measurements, 
seeds were surface sterilized and stratified at 4 °C for 3 d. Seeds were 
sown on 80% Murashige and Skoog (MS) medium plus vitamins, 
1% sucrose and kanamycin (35mg I" 1 ) and placed in a growth cham- 
ber (ATC 26, Controlled Environments; photosynthetic photon flux 
(PPF) of approximately 120-150 umol nr 2 s" 1 ; 10 h:14h lightdark; 
22 °C: 19 °C day:night). After 6-10 d of growth on selection medium, 
seedlings were transplanted into autoclaved Promix soil and grown 
in an ATC 26growth chamber as above. After 1 week of growth on 
soil, plants were given weekly fertilizer treatments with 0.4 or 0.8 g 
H of Peter's fertilizer. Plants were about 30^12 d old when used for 
experiments. 

To examine the endogenous expression of ATHB17 transcripts 
by in situ RNA-RNA hybridization, seeds were prepared as above 
and germinated on solid 80% MS medium (80% MS plus vitamins, 
0.3% sucrose, and 1% BactoAgar) with a 16 h:8h light:dark pho- 
toperiod for 6 d. Seedlings were then transferred to vertical plates 
and grown for an additional 6 d. Seedlings that were 12 d old were 
used for in situ analysis of shoot apices and roots. Soil-grown 
plants were used to analyse the later stages of Arabidopsis growth. 
Seeds were germinated directly in soil after 2 d of stratification 
in the dark at 4 °C. Plants were grown in a 16 h:8h light:dark 
photoperiod at 22 °C in a growth chamber. Inflorescence, floral 
meristems, and developing seeds were collected from 4-5-week- 
old plants. 

For laser-capture microdissection, 11-d-old seedlings (grown as 
above for in situ analysis) from athbl7-l and its wild-type sibling 
were used for collection of specific root tissues. 



RNA isolation and quantitative RT-PCR (qRT-PCR) 

Expression of the ATHB17 mRNA transcripts in athbl7 mutant 
and overexpressing plants was examined by qRT-PCR. RNA 
was isolated from frozen root tips or 14-day-old seedlings using a 
Qiagen RNeasy Plant Mini kit and reverse transcribed (Superscript 
II RT; Invitrogen) using a mixture of oligo(dT) and random hex- 
amer primers (60 and 40% respectively). cDNA templates were 
then analysed by real-time PCR (7900 HT Fast Real-Time PCR; 
Applied Biosystems) using a SYBR Green PCR Mix (Applied 
Biosystems). 

Primers used for the RT-PCR were: ATHB1 7 transformation con- 
struct (S'-TTCGTGTAAATACTAAGAGACTCTGTTCCG-S' and 
5'-TGCCATAATACTCAAACTCAGTAGGA-3'), ATHB17 native 
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expression (5'-AGGTGGTTTGGTTCATTAACGGAAG-3' and 
5'-GCAGCGAGGACACATAGTAAGGC-3'),^r//52(5'-ACATG 
AGCCCACCCACTACTTTGAC-3' and 5'-CAGGAGCCCACGC 
ATTGACC-3'),^7W^72(5'-AGCACAATACTCTCAATCCCA-3' 
and 5'-CTGCTTTAACTTTGTCCTTGCTC3'), A THA Tl (5'-GAA 
CACAACACTCTCAATCCC-3' and 5'-TTCAAGTATTCGCAAT 
CTACCTCC-3') and a mid-level constitutive control gene AT2G32170 
(5'-GCAGGAATAACAGAAGGTTTCTCTATGTGT-3'and5'-AT 
GAGCAGTGTCGATGAAAAAGCA-3'; Czechowski et al, 2005). 
Expression levels were calculated based on AC t , the difference in 
cycle threshold (C t ) between the gene of interest and the control, 
with signals normalized to AT2G32170, a mid-level constitutive 
control gene. 

Protoplast-based transcriptional repression assay 

The cDNA sequence encoding the ATHB17 protein was cloned 
downstream of the cauliflower mosaic virus 35S promoter and a 
translational enhancer from the 5'-leader of tobacco mosaic virus 
(Skuzeski et al, 1990). All promoter sequences were PCR amplified 
and cloned upstream of the GUS reporter gene. The Agrobacterium 
NOS transcriptional terminator was used in all reporter and effec- 
tor plasmids (Skuzeski et al, 1990). The sequence length of each 
promoter was as follows: ATHAT1 (1821 bp), ATHB2 (1500 bp), 
and AGL8 (1833bp). The nucleotide position of each promoter was 
the number of bases upstream of the translation start site of the 
respective genes. The 35S::CAT reporter gene construct encoding 
chloramphenicol acetyltransferase (CAT) has been described previ- 
ously (Tiwari et al. , 2006). 

Protoplasts were prepared from leaves of 3^1-week-old 
Arabidopsis plants using cellulase (Research Products International) 
and macerozyme (SERVA Electrophoresis; Tiwari et al., 2006). 
For the assay, 200 ul of protoplast cells was co-transfected with 
10 ng °f the ATHATlv.GUS, ATHB2::GUS, or AGL8..GUS 
reporter gene construct and either 10 |xg of the 35 S.. CAT control 
construct or a 35S:.ATHB17 construct. Plasmid DNA (EndoFree 
Plasmid Maxi; Qiagen) was introduced into protoplasts by incu- 
bating with 40% polyethylene glycol for 30min. After polyethylene 
glycol removal, protoplasts were kept in darkness at room tem- 
perature for 18-20h. Cells were lysed and incubated with ImM 
4-methylumbelliferyl-(3-D-glucuronide (Gold Biotechnology) at 
37 °C for 1 h. P-Glucuronidase (GUS) activity was quantified by 
taking fluorescence measurements using a Synergy HT Microplate 
Reader (BioTek). GUS measurements for each construct were 
averaged (rc=3). 



In situ RNA-RNA hybridization 

DNA templates were generated by PCR using gene-specific primers 
(A THB1 7-forward, 5'-CTTCGATCCTAGCTCTTAAGAACC-3', 
and^ri/fi77-reverse,5'-GTAACCGATTCATGTCTAGCC-3')with 
incorporation of a T7 promoter sequence (5'-CTAATACGACTCA 
CTATAGGG-3') on one primer (Baklanov et al, 1996; Iacobuzio- 
Donahue et al, 2002). To generate template for the sense and anti- 
sense probes, the T7 promoter sequence was incorporated at the 
5' end of A THB1 7-forward and the 5' end of ATHBl 7- reverse, 
respectively. 

Digoxigenin (DIG)-labelled RNA probes were generated by in 
vitro transcription using T7 RNA polymerase (Life Technologies) 
except that a 1:2 ratio of UTP:DIG-11-UTP (Roche Applied 
Science) was used during the reaction. Fixation of tissues, prepa- 
ration of sections, hybridization, and washes were carried out as 
described by Jeffrey Long (http://www.its.caltech.edu/~plantlab/ 
protocols/insitu.pdf, last accessed 26 July 2013) with the following 
modifications: (i) RNAase treatment was omitted; (ii) the blocking 
reagent was 1% blocking reagent in lx maleic acid buffer (Roche 
Applied Science); (iii) the wash buffer contained lx maleic acid 
buffer with 0.03-0.05% Tween 20 (v/v) (Roche Applied Science); and 
(iv) sections were mounted using Immu-Mount (Thermo Scientific). 



Preparation of paraffin sections and laser-capture 
microdissection 

Primary and lateral root tips from 11-d-old seedlings were used 
for the collection of specific root tissues. Preparation of paraffin- 
embedded tissues (Kerk et al, 2003; Jiao et al, 2009) was per- 
formed with minor modifications. Specifically, root tips were fixed 
in ice-cold fixative (70% ethanol, 15% acetic acid, 0.6% trehalose) 
for approximately 3 h, and then treated twice with 75% alcohol for 
30min each and 90% alcohol for 45min. Stepwise infiltration of 
paraffin was done with a TP1020 automatic tissue processor (Leica 
Microsystems): 100% ethanol for lh, 0.2% Eosin Y (in 100% eth- 
anol) for 30min, 100% ethanol for 2h, ethanol/Histo-Clear series 
(3:1, 1:1, 1:3, 0:1, 0:1, 1-2 h each), and paraffin (for 4 and 8 h with 
vacuum). Paraffin sections (8 um) were prepared as described previ- 
ously (Jiao et al, 2009). 

Paraffin sections were deparaffinized in Histo-Clear and air 
dried before slides were viewed using the fluorescence mode and 
the rhodamine reflector of a PALM MicroBeam System (Carl 
Zeiss) to facilitate the identification of different cell types. The 
root cap, quiescent centre (QC), epidermis/endodermis/cortex, and 
stele tissues were collected into adhesive caps (Carl Zeiss). Three 
biological replicates, each pooled from 40-60 roots, were collected 
for each sample 

Microarray analysis 

Total RNA was prepared and amplified and cRNA was gen- 
erated from microdissected tissues using commercial kits 
(PicroPure RNA, RiboAmp HS RNA, and Turbo Labeling Biotin 
kits, respectively; Life Technologies). Transcript profiling experi- 
ments were performed using a custom full-genome Arabidopsis 
Affymetrix GeneChip (mbiAThla520184) microarray and ana- 
lysed according to Holtan et al (2011). Calculations for the over- 
representation analysis were performed according to Holtan et al 
(2011). 

The Affymetrix data files (.CEL), pre-processed profiles, and 
supplementary files described herein have been deposited at the 
National Center for Biotechnology Information (http://www.ncbi. 
nlm.nih.gov/geo, last accessed 26 July 2013) under GEO accession 
no. GSE48127. 



Preparation and analysis of mesophyll cell suspensions for 
chloroplast counting 

Chloroplasts were counted using Nomarski differential interference 
contrast optics (Leica DMR-HC; Leica Microsystems) in individual 
fixed mesophyll cells obtained by the maceration of prepared leaf 
tissue (Pyke and Leech, 1991). 



Estimating biochemical limitations to light-saturated 
photosynthesis 

The protocol of Long and Bernacchi (2003) guided the gen- 
eration of plots of light-saturated photosynthesis (A sllt ) and the 
rate of linear electron transport through photosystem II (/ PSII ) 
to substomatal [C0 2 ] (Q) for the youngest, fully expanded 
leaves of plants for a control line and two 35S..ATHB17 lines 
using an infra-red gas analyser with a modulated chlorophyll 
a fluorimeter (LI-6400 XT; Licor). Initial measurements were 
made at an atmospheric [C0 2 ] (C a ) of 400 (xmol of C0 2 moL 1 , 
after 40min exposure to a PPF of 700 u.mol m 2 s _1 . After the 
initial measurement, further measurements were made after 
three stepwise decreases in C a to 150, 100, and 50 umol of C0 2 
moL 1 . After a second measurement at 400 umol of C0 2 moL 1 , 
three final measurements were made after stepwise increases in 
C a to 600, 700, and 800 umol of C0 2 moL 1 . Measurements for all 
three lines were made at the same leaf-air water vapour pressure 
deficit of about 1 kPa and relative humidity of about 50%. 
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The curve-fitting software of Sharkey et al. (2007) and kinetic 
constants described therein were used to make in vivo estimates 
of: (i) the maximum rate of ribulose-l,5-bisphosphate (RuBP)- 
saturated photosynthesis ( V c max ), from the three measurements of 
A sal and C; made at the lowest C a levels; and (ii) the light-saturated 
capacity for RuBP regeneration, calculated and expressed in terms 
of electron flow required to support the modelled rate of A sM (/ max ), 
from the three highest C a levels. Estimates of / PSII were estimated 
as the product of the operating efficiency of PSII and the fraction 
of incident PAR absorbed by PSII (Genty et al. 1989). The frac- 
tion of the 700 umol rar 2 s" 1 of PAR incident upon the leaf in the 
leaf chamber that was transmitted through leaves was measured 
by pressing a quantum sensor (LI 190; Licor) to the centre of the 
abaxial side of the leaf. No measure of reflectance from the adaxial 
leaf surface was made; consequently absorption was calculated as 
1 - transmission for the determination of / PS n- We assumed that 
50% of incident PAR was absorbed by PSII for the purposes of 
calculating / PSII . 



Chlorophyll assays 

Chlorophyll content was measured using a SPAD 502 (Konica 
Minolta). 

Data analysis 

Mesophyll cell plan areas, chloroplast number per cell, and indi- 
vidual chloroplast plan areas were measured directly from the 
microscope image using ImageJ (Abramoff et al, 2004). Data 
from leaves were averaged to give a mean and standard error (SE). 
Significance was tested with a Student's /-test and a standard lin- 
ear regression model using R (version 2.2.0; R Development Core 
Team, 2007). 

The hypothesis that 35S..ATHB17 plants had increased K Cjmax 
and / max compared with the empty-vector control plants was tested 
with a one-tailed Student's /-test using Excel (Microsoft). Data 
analysis was performed on six replicate samples unless otherwise 
stated and an effect was described as statistically significant when 
P< 0.05. 



Results and discussion 

ATHB1 7 is expressed in the QC of the primary and 
lateral roots 

We analysed the spatial expression pattern of ATHB17 by 
in situ RNA-RNA hybridization (Fig. 1A-F) using DIG- 
labelled antisense and sense probes targeted to a unique 
region of the ATHB17 coding sequence (nt 173-353). Of 
the tissues examined (roots, shoot apical meristems, inflores- 
cence meristems, floral meristems, and developing embryos), 
hybridization was only detected in specific cells in embry- 
onic, primary, and lateral roots. In a survey across a devel- 
opmental series, ATHB17 signal was first detected in the 
base of globular stage embryos, and was then restricted at 
the transition stage to the lens-shaped cell, which gives rise 
to the QC. In serial sections from older embryos and the 
primary root tip of seedlings (Fig. 1A-E) as well as the tips 
of lateral roots (Fig. IF), ATHB1 7 expression continued to 
be localized specifically to the small number of cells that 
comprised the QC. No signal was detected using the control 
sense probe (Fig. 1G-K). 

The root is comprised of several tissues that undergo pre- 
dictable divisions. Initials for all root tissues surround a group 
of four to seven mitotically less active cells known as the QC. 




Fig. 1. In situ hybridization using DIG-labelled antisense and 
sense probe targeting nt 1 73-353 of the ATHB17 coding 
sequence. (A-F) Antisense probe, (G-K) sense probe. (A-D, G-J) 
Embryos at the globular stage (A, G), transition stage (B, H), 
heart stage (C, I), and linear cotyledon stage (D, J). (E, K) Primary 
root tip. (F) Lateral root tip. In each panel, the region of greatest 
hybridization signal is highlighted with a black arrow. Bars, 
50 urn. 



The QC is important for proper root development and regu- 
lates processes and events in neighbouring cells for proper 
meristem establishment and maintenance (Doerner, 1998; 
Jiang and Feldman, 2005). The root QC-specific expression 
of ATHB17 implied that ATHB17 plays a role in the regula- 
tion of root growth and development. 
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Comparison of root QC transcript profiles between 
athbl 7-1 and its wild-type sibling 

We screened T-DNA insertion lines (Alonso et al, 2003) and 
identified three different ATHB17 alleles (Supplementary 
Fig. SI at JXB online). RT-PCR analysis indicated a signifi- 
cant reduction of the normal ATHB17 transcript in the line 
designated athbl 7-1, containing a T-DNA insertion of 162 bp 
into the third intron of ATHB17 (Supplementary Fig. S2A at 
JXB online). Western blot analysis confirmed the absence of 
ATHB17 in the root tips of athbl 7-1, but not in athbl 7-2 and 
athbl 7-3 (Supplementary Fig. S2B). The athbl 7-1 line was 
used for further analysis. 

We examined how ATHB17 regulates gene expression by 
transcript profiling of athbl 7-1 and its wild-type sibling. 
To increase our sensitivity to effects caused by the loss of 
ATHB17, we analysed tissue samples (collected by laser- 
capture microdissection) that were highly enriched with QC 
cells. To validate that we were able to enrich for QC cells by 
laser microdissection, a small panel of genes was selected 
and their transcript levels analysed by qRT-PCR across four 
distinct tissue types dissected from the root tip: QC, root cap, 
stele, and epidermis. The expression of known QC-localized 
genes [PLETHORA 1 (PLT1), WIP4, and ATMYB34; Aida 
et al., 2004; Nawy et al., 2005] was enriched in QC sam- 
ples (Supplementary Fig. S3A at JXB online). On the other 
hand, genes known to be expressed in stele (AT5G25490; 
Nawy et al, 2005) or root cap tissue (ATMDK-20; del 
Campillo et al, 2004; Lilley et al, 2011) displayed stele- or 
root cap-enriched expression, respectively (Supplementary 
Fig. S3B, C). 

Of the three class II HD-Zip genes expressed in the root 
QC (Nawy et al, 2005), ATHB2 and ATHAT2 showed 
increased expression in the QC of athbl 7-1 roots relative 
to the wild-type control, while ATHAT1 did not (Fig. 2A). 
We confirmed the effect on ATHB2 and lack of effect on 
AT H ATI by performing qRT-PCR using RNA isolated 
from whole root tips of athbl 7-1 and its wild-type sibling 
(Fig. 2B). The more subtle effect on ATHAT2 observed in 
the array experiments could not be confirmed in whole root 
tips. The array and qRT-PCR results suggested that ATHB17 
functions to repress the expression of ATHB2 and ATHAT2 
in the QC, or that ATHB17 may act to repress an activa- 
tor of these genes. Indeed, in seedlings of 35S::ATHB17 
lines in which ATHB17 was constitutively expressed at 
high levels, ATHB2, ATHAT1, and ATHAT2 were all 
strongly repressed (Fig. 2C). We tested whether ATHB17 
could regulate the expression of ATHB2 and ATHAT1 
using a transient expression assay. We co-transfected con- 
structs encoding ATHB17 or CAT and pATHB2::GUS or 
pATHATT.GUS reporter genes (Fig. 2D) into Arabidopsis 
mesophyll protoplasts. The pAGL8::GUS reporter gene 
was used as a non-specific reporter that was not expected 
to be regulated by ATHB17 protein. The transcription of 
the reporter gene driven by either class II HD-Zip pro- 
moter was strongly repressed in the presence of ATHB17. 
Thus, class II HD-Zip gene expression within the root QC 
was regulated by other class II HD-Zip proteins and may 



function as a negative regulatory loop in which homeostasis 
is mutually controlled (Ariel et al, 2007; Sorin et al, 2009; 
Harris et al, 2011). The specific expression of ATHB17 in 
the root QC suggested that ATHB17 may be involved in 
some aspects of root development. However, the roles that 
ATHB2, ATHAT1, ATHAT2, and ATHB17 play in affect- 
ing root growth are unknown at this time. 

Several related biological processes are subtly 
perturbed in athbl 7-1 

Gene set over-representation analysis using gene ontology 
(GO) terms was used to determine which biological pro- 
cesses and molecular functions showed the greatest differ- 
ences between athbl 7-1 and its wild-type sibling (Fig. 3). 
Genes differentially expressed between athbl 7-1 and its 
wild-type sibling with annotation to biological processes 
that included responses to auxin stimulus, embryo develop- 
ment, and axis polarity were over-represented. Functional 
analyses indicate that some HD-Zip class II genes (ATHB2, 
ATHB4, and ATHAT2) are linked to auxin (Steindler et al, 
1999; Sawa et al, 2002), consistent with our observations 
on the expression of class II HD-Zip genes in athbl 7-1 
QC cells. 

The presence of genes in GO categories related to embryo 
or seed development and axis polarity (Fig. 3) was interest- 
ing considering that, while class II HD-Zips are implicated 
in the integration of shade-avoidance responses and hor- 
mone-mediated growth (Sawa et al, 2002; Sorin et al, 2009), 
class III HD-Zip proteins and KANADI are characterized 
as developmental regulators during embryogenesis, the estab- 
lishment of the apical meristem, lateral root formation, and 
the determination of leaf polarity (Hawker and Bowman, 
2004; Prigge et al, 2005; Izhaki and Bowman, 2007). The 
expression of PHAVOLUTA (AT1G30490), KANADI 2 
(AT1G32240), and REVOLUTA (AT5G60690) was induced 
in the QC of athbl 7-1 roots (Supplementary Table S2 at JXB 
online). A common link between class II/III HD-Zip proteins 
and KANADI-related proteins is auxin. Auxin concentration 
and flux mediated by PIN-FORMED1 (PIN1) contribute to 
embryo patterning during embryogenesis and the establish- 
ment of the apical-basal axis and bilateral symmetry (Izhaki 
and Bowman, 2007). 

A number of genes with GO terms specifying plastid- 
related processes and constituent parts including thylakoid 
membrane organization, stroma, stromule (Kohler and 
Hanson, 2000), and the term 'plastid part' were significantly 
over-represented, with more genes repressed than induced in 
the athbl 7-1 mutant. 

A whole-plant phenotype comparison among wild-type, 
35S::ATHB17, and athbl 7-1 Arabidopsis plants 

The athbl 7-1 plants displayed no clear differences from 
wild-type controls when grown on plates to examine poten- 
tial root phenotypes or when grown in soil to examine over- 
all morphology and chlorophyll content (Supplementary 
Fig. S4 at JXB online), which suggested that ATHB17 is 
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Fig. 2. Expression levels oiATHAT2, ATHAT1 , and ATHB2 in various Arabidopsis tissues. Microarray analysis indicated that the 
transcript levels of selected class II HB-Zip genes were altered in the QC-enriched samples ofathb17-1 compared with its wild-type 
sibling. (A) ATHAT2 and ATHB2 expression, but not ATHAT1 expression, was induced in the athb17-1 QC. (B) qRT-PCR was used 
to verify the altered expression pattern observed in panel A. (C) ATHAT1 , ATHAT2, and ATHB2 expression is significantly repressed 
in seedlings of two 35SvATHB17 lines (line 361 and line 370). (D) ATHB1 7 also repressed the transcription ofATHATI and ATHB2 in 
Arabidopsis mesophyll protoplasts. Full-length ATHB 1 7 or CAT was co-transfected with the pATHATl : :GUS or pATHB2::GUS reporter 
gene into Arabidopsis mesophyll protoplasts. The pAGL8::GUS reporter gene was used as a non-specific reporter as it was not 
expected to be regulated by ATHB1 7. The data are means ±standard deviation (SD) of three replicates. PCR comparisons were made 
after normalization to a constitutive control gene by the ACt method. 



functionally redundant with other class II HD-Zip proteins 
in the root QC. The primary and branch root lengths as well 
as branch root number were similar in both athbl7-l and its 
wild-type sibling (Supplementary Fig. S4A-C). No detect- 
able rosette phenotype was observed (Supplementary Fig. 
S4D, E). On the other hand, athbl7-l may produce a con- 
ditional phenotype or our experimental growth conditions 
may not have been sensitive enough to detect a discernible 
phenotype. 

In contrast, overexpression of ATHB1 7 produced marked 
phenotypic effects; in particular, the leaves of the transgenic 
lines were consistently visibly darker green (Fig. 4A). A series 
of 35S::ATHB17 lines was examined for this phenotype; 



chlorophyll content as measured by soil plant analysis devel- 
opment (SPAD) was higher in all lines (Fig. 4B) and was pos- 
itively correlated with ATHB17 mRNA levels. Chlorophyll 
levels in leaves as measured by SPAD from athbl7-l plants 
were comparable to those seen in empty-vector control 
leaves (Supplementary Fig. S4F). Leaf blade development in 
35S::ATHB17 lines was also de-repressed along the petiole, 
creating a loss of distinct proximal and distal zones and result- 
ing in a slight leaf blade on petiole phenotype that resulted in 
a rounder appearance (Fig. 4A; Ha et al, 2003; Hepworth 
et al, 2005). We have suggested that, within the root QC, 
ATHB 17 with ATHB2, ATHAT1, and ATHAT2 form a neg- 
ative regulatory loop It has been widely demonstrated that 
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Fig. 3. Functional category analysis of genes differentially expressed in the root QC ofathb17-1 compared with the wild-type sibling 
line. Over-representation of terms from different aspects of GO, including (A) Biological Process, (B) Molecular Function, and (C) Cellular 
Compartment are shown separately. Based on analysis of all GO terms annotated to Arabidopsis genes (Berardini ef al. , 2004), the 
categories listed above were over-represented [P <0.01) among significantly mis-expressed genes. For each over-represented GO 
category, we have shown the number of genes with significantly altered expression and the total number of genes annotated to the 
term. The horizontal length of the bar indicates the ratio of the number of genes annotated to the term that have significant expression 
changes relative to the number of genes from that category that would have been expected to show such changes by chance, while the 
vertical length of the bar represents the magnitude of the average absolute fold change of these genes. The percentage of induced and 
repressed genes annotated to each GO term is indicated by red and green, respectively. See Supplementary Table S1 at JXB online for 
a list of genes over-represented in each of the GO categories. 



the expression of class II HD-Zip genes is mutually regulated, 
in that when one class II HD-ZIP gene is overexpressed, the 
expression of other class II HD-Zip genes is repressed (Ohgishi 
et al, 2001; Sawa et al, 2002; Sorin et al, 2009; Harris et al, 
201 1). The phenotypes observed with constitutive overexpres- 
sion of ATHB1 7 could result from repression of other class II 
HD-Zip genes or other native targets during embryo develop- 
ment and growth during the seedling and rosette stages. The 
slight blade on petiole phenotype (Ha et al, 2003; Hepworth 
et al, 2005) suggested that subtle alterations in leaf devel- 
opment are caused by ATHB17 overexpression. Leaf pattern 
formation is dependent on the establishment of the abaxial/ 
adaxial axis because loss of either abaxial or adaxial identity 
perturbs leaf blade outgrowth and symmetry. Consequently, 
some leaf phenotypes in 35S::ATHB17 plants may result 
from perturbations in abaxial/adaxial polarity mediated by 
effects on auxin transport/signalling or the disruption of a 
complex regulatory network consisting of various class II 
and class III HD-Zips plus KANADI-like proteins (Harris 
et al, 2011; Brandt et al, 2012). 

Chloroplast and cell characteristics of 
35S::ATHB1 7 plants 

Representative cellular phenotypes of typical mature leaf 
mesophyll cells of an empty-vector control line and three 



35S::ATHB17 lines (306, 314, and 319; T2 generation) 
are illustrated in Fig. 5A. Control and 3 5 S:\ATHB17- 
overexpressing lines showed a close relationship between 
plastid number and mesophyll cell size (Fig. 5B). 
35S::ATHB17 plants accumulated more plastids per unit 
cell area compared with the control line (Table 1). The dis- 
tribution of chloroplast plan area in fully expanded mes- 
ophyll cells was also investigated and is shown in Table 1. 
35S::ATHB17 lines exhibited a slight reduction in plastid 
plan area compared with the control line. The proportion 
of mesophyll surface area covered by chloroplasts was con- 
sistent between 55S::;477/,Bi7-overexpressing lines and the 
control line (data not shown). This may result from compen- 
sation between chloroplast number and chloroplast size; i.e. 
35S::ATHB17 mesophyll cells have more chloroplasts per 
cell than a wild-type control line, but the chloroplasts are 
somewhat smaller. 

Based on published data (Pyke and Leech, 1992; Pyke, 
1997), the smallest mesophyll cells (500-800 um 2 , representa- 
tive of post-meristematic cells) of fully expanded mature 
leaves have approximately 10-20 chloroplasts. However, based 
on our regression analysis (Fig. 5B), the number of plastids 
in the smallest cells of 35S::ATHB1 7 plants was greater than 
that found in the youngest cells of an empty-vector control 
line. Therefore, the number of proplastids per meristematic 
cell is potentially greater in 35S::ATHB17 plants than in 
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Fig. 4. (A) Plants (30 d old) grown under short-day conditions illustrating that the transgenic plants (T3) had leaves that were darker 
green, and rounder with a slight leaf blade on petiole phenotype.(B) A high correlation was found between chlorophyll content 
determined by SPAD (blue) and ATHB17 transgene expression (black, shown relative to the constitutive control gene AT2G32170) 
was demonstrated across 1 1 T2 35S::ATHB1 7 lines. Results are shown as the mean±SE (n=20). **P <0.05. 



controls (Table 1). These data suggest that ATHB17 overex- 
pression regulates proplastid division rather than chloroplast 
division. 

Light-saturated photosynthetic capacity is increased in 
35S::ATHB1 7 Arabidopsis plants 

Two well-characterized 3 5 S::ATHB1 7 lines (373-3 and 378-1; 
T3 generation) were studied to investigate whether increases 
in leaf chlorophyll content would impact photosynthesis in 
Arabidopsis. For both 35S::ATHB17 lines, leaf SPAD was 
significantly increased, by 36% relative to the control line 
(Table 2). Increased leaf chlorophyll was consistent with a 
significant 6% decrease in transmission of PAR through the 
leaves in the 35S::ATHB17 lines (Table 2). The decrease in 
transmission of PAR through the 35S::ATHB17 leaves was 
probably due to increased PAR absorption, an effect that 
could be reconciled with a significant increase in / PSII of up 
to 14% in the two 35S;:ATHB17 lines, relative to the control 
line, over a range of Q (Fig. 6A). This increase in / PSII was 
largely the consequence of increased light absorption, with 
increases in the operating efficiency of PSII (F q /F m ) in the 



35 S:\ATHB17 lines typically being smaller over the same 
range of Q (Fig. 6A). 

The plots of /p SII and A mt (Fig. 6B) against increasing 
C- provided evidence of increased photosynthetic capac- 
ity in both 35S::ATHB17 lines. A progressive stimulation 
of y4 sat relative to the control line was seen as Q increased 
in both 35S::ATHB17 lines (Fig. 6B). AT the three high- 
est Q values, both A S!it and / PSII were significantly higher 
in the 35S::ATHB17 lines; however, while A S&1 was still ris- 
ing in all lines, / PSII was reaching, or had reached, a plateau. 
Increasing A sM over a range of Q where / PSII is at a plateau 
is indicative of A sM being limited by the capacity to regener- 
ate RuBP in the Calvin cycle (Long and Bernacchi, 2003). 
For both 35S::ATHB17 lines, / max estimated from these 
three data points was significantly higher than for the con- 
trol line (Table 2). At the three lowest Q values, A sat and J PSII 
increased linearly with increasing Q, observations indicative 
of A S3l being limited by the activity of ribulose-l,5-bispho- 
sphate carboxylase oxygenase (Rubisco; F cmax ). For both 
35S: :A THB1 7 lines, F c max estimated from these three points 
was also significantly higher than in the control line (Fig. 6 
and Table 2). Stomatal conductance, measured at ambient C a , 
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Fig. 5. Isolated leaf mesophyll cells (A) and chloroplast number per mesophyll cell versus cell size (B) from the youngest fully expanded 
leaves of 40-d-old control or 35SvATHB17 lines. (A) Mesophyll cells were viewed with Nomarski differential interference contrast optics. Bar, 
25 (im. (B) Each data point represents the measurement from one cell. The solid red line in each figure represents the regression line for the 
empty-vector control line, while the solid black line is the regression of the transgenic in each panel. Values for r 2 are indicated on the graphs. 



Table 1 . Mean chloroplast plan area, mean mesophyll plan area, and plastid number for populations of mesophyll cells from the 
youngest fully expanded leaves of 40-d-old control and 35S::ATHB1 7 Arabidopsis plants 

Chloroplast plan areas were measured by image analysis of fixed isolated cells from three 35S::ATHB 1 7-overexpressing lines (T2 generation) 
and an empty-vector control line (pMEN65) in two independent experiments. Plastid number per cell (±SE) was determined from a regression 
line of chloroplast number per cell on mesophyll cell plan area using a value of 1 000 |im 2 for the mesophyll cell plan area (Fig. 5). The y-intercept 
of the regression line was used to estimate the number of proplastids (+SE) in very small, meristematic cells. For Experiments 1 and 2, mean 
chloroplast plan area is the mean +SE of 600 chloroplasts from 30 or 60 different mesophyll cells, respectively. **P <0.05. 



Genotype 


Mean chloroplast 
plan area (|im 2 ) 


Plastids 




No. chloroplasts per 1000 pm 2 
mesophyll cell plan area 


Estimated number 

proplastids per meristematic cell 


Experiment 1 








pMEN65 


39.8 + 0.4 


35 + 2 


14 + 2 


35S::ATHB17 line 306 


35.9 + 0.3** 


51 + 1** 


33 + 1** 


35S-.-ATHB17 line 314 


35.0 + 0.3** 


55+1** 


37 + 1** 


35S-.-ATHB17 line 319 


38.3 + 0.4** 


48 + 1** 


25 + 1** 


Experiment 2 








pMEN65 


36.7 + 0.4 


44+1 


22 + 2 


35S::ATHB1 7 line 306 


32.7 + 0.3** 


70+1** 


48 + 1** 


35S::ATHB17 line 314 


32.7 + 0.3** 


67+1** 


44 + 1** 


35S-.-ATHB17 line 319 


30.5 + 0.4** 


57+1** 


34 + 1** 



Table 2. Key components of photosynthetic capacity are increased in 35S::ATHB1 7 lines 

Sample means ±SE are shown for leaf SPAD, percentage transmission of PAR, and estimates of V crrm and J wax . Data were collected on at least five 
replicate plants for the empty-vector control and two 35S:ATHB1 7 lines. The means of each data set collected for the two 35S:ATHB1 7 lines were 
significantly different from the control line (P <0.05; two-tailed Student's f-test). 



Genotype SPAD Transmission (%) V CJIlax J max 

pMEN65 31.2 + 0.5 12.7 + 0.5 28.8 + 1.1 80.4 + 2.3 

35S::ATHB1 7 line 373-3 43.9 + 2.0 5.4 + 0.5 33.1 + 1.6 90.2 + 2.3 

35S-.-ATHB17 line 378-1 42.4 + 1.4 5.8 + 0.4 33.6 + 2.1 87.5 + 2.7 
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Fig. 6. Plot of Jpsn and F v /F m (A) and A sat (B, C) against Q for 
select 38-d-old 35S::ATHB17 lines and an empty-vector control 
line (pMEN65). The data shown in (C) give an expanded view of 
the C0 2 -Iimited part of the larger A sat /Q response curve shown 
in (B). The photosynthetic phenotypes reported have proved 
to be highly penetrant across multiple 35S::ATHB 1 7 lines, and 
the effects described were repeatable in multiple independent 
experiments. Results are shown as the mean±SE (n=6) from a 
representative experiment. 



in the Calvin cycle, and conductance of C0 2 to the sites of 
carboxylation are all key components of the overall strat- 
egy (Long et al., 2006; Zhu et al, 2010; Raines, 2011). This 
study provides evidence that modulation of A THB17 activity 
can increase the biochemical capacity for photosynthesis in 
Arabidopsis, and could improve leaf photosynthesis at cur- 
rent and predicted future atmospheric C0 2 concentrations. 

Conclusion 

In the QC, ATHB17 appears to be part of a complex negative 
regulatory loop with other class II HD-Zip genes. Based on 
microarray data, this regulatory loop controls processes such 
as auxin transport, embryo development, and axis polarity. 
Constitutive overexpression of ATHB17, on the other hand, 
modulates leaf development, producing changes in leaf mor- 
phology, plastid number and size, and the capacity for pho- 
tosynthesis. Such results underscore the utility of functional 
genomics research to identify candidates for the modification 
of complex traits in plants. 



Supplementary data 

Supplementary data are available at JXB online. 

Supplementary Fig. SI. Loss-of-function mutants for 
ATHB17. 

Supplementary Fig. S2. Expression levels of ATHB17 
in the corresponding T-DNA insertion lines represented in 
Supplementary Fig. SI. 

Supplementary Fig. S3. Quantitative RT-PCR analysis 
to analyse the differential expression of various genes using 
RNA isolated from the root tip, quiescent centre (QC), 
and stele. 

Supplementary Fig. S4. When grown under our experi- 
mental conditions, athbl7-l does not produce a visible 
phenotype. 

Supplementary Table SI. Genes found to be significantly 
regulated, where the Benjamini-Hochberg-corrected P value 
is less than 0.05 and the fold change is greater than 1.3, 
are shown. 

Supplementary Table S2. Genes found to be significantly 
regulated, which underlie the over-representation presented 
in Fig. 3, are shown. 
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was reduced by 17% in the 35S::ATHB17 373 line and by 9% 
in the 35S::ATHB17 378 line; however, neither effect was sta- 
tistically significant (P >0.2; data not shown). 

Much research effort is currently focused on increasing 
the C3 photosynthetic capacity. Increasing the maximum 
activation state of Rubisco, the capacity to regenerate RuBP 
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